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The three-phase contact line formed by the intersection of a liquid-vapor interface of an electrolyte 
solution with a charged planar substrate is studied in terms of classical density functional theory 
applied to a lattice model. The influence of the substrate charge density and of the ionic strength 
of the solution on the intrinsic structure of the three-phase contact line and on the corresponding 
line tension is analyzed. We find a negative line tension for all values of the surface charge density 
and of the ionic strength considered. The strength of the line tension decreases upon decreasing the 
contact angle via varying either the temperature or the substrate charge density. 
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I. INTRODUCTION 

The line tension is the free energy per length associ¬ 
ated with the contact line where three phases meet in 
space. For example, for a sessile liquid drop on a solid 
substrate and surrounded by the vapor phase, the con¬ 
tact line corresponds to the periphery of the circle where 
the liquid-vapor interface meets the substrate. Although 
its magnitude is small (both theoretical and emerimental 
values are of the order of 10“^^ to 10“^^ N the line 

tension plays an important role for various^stems and 
phenomena such as spreadin g o f droplets [l0|, UJj, wet¬ 
ting of nanoporons surface s stability of emulsions 
and foams [l^, drop size [ij, and many others. The 
line tension has been the subject of numerous theoretical 
and experimental investigations (see, e.g.. Refs. [iEHii 
and references therein). Experimental setups to study 
line tensions encompass solid- l^uid-gas systems, such as 
drops on solid substrates [El, [Ef, bubbles on solid sub¬ 
strates or on particles at liquid-gas interfaces [ni- 

[ 2 ^, and liquid-liquid- va^r systems, such as liquid lenses 
at liquid-gas interfaces ]l|, 0, [13 • Theoretical investiga¬ 
tions include extensions of capillarit y th eory pEj . which 
take into account line tension effects [^, , microscopic 

theories [H, , as well as molecular dynamics [ol, Isl - 

[s?! and Monte Carlo [1, [H, [39| simulations. Except for 
simulations of pure water, these inyestigations deal with 
simple fluids or binary liquid mixtures thereof. Howeyer, 
most actual fluids used in wetting applications comprise 
seyeral components. Eor polar fluids, such as water, these 
additional components often carry an ionic character. It 
is well-known that the presence of ions in a fluid creates 
the Debye length as an additional length scale, which in¬ 
creases upon decreasing the ionic strength. Eor dilute 
electrolyte solutions it typically exceeds the bulk corre¬ 
lation length of the pure solyent. Hence the structure 
and the wetting behayior of an electrolyte solution can 
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be expected to differ significantly from that of the pure 
solyent jiol-lill. 

There are only few studies concerned with the influence 
of electrostatic interactions on the line tension jiEI-Esj. 
In Ref. [13 the theory of capillarity has been extended 
taking into account line contributions as well as electric 
charges at the interfaces and at the three-phase contact 
line (TPCL). Within this approach electrowetting has 
been interpreted as a line tension effect, but some of the 
corresponding predictions are in disagreement with ex¬ 
perimental data dil . In Ref. Q an equation for the con¬ 
tact angle as function of the electrostatic potential at the 
TPCL and an estimate for the electrostatic contribution 
to the line tension haye been deriyed using a yariational 
approach for a wedge-like geometry. Based on a Poisson- 
Boltzmann theory the analysis in Refs. 111113 considers 
only the electrostatic part of the free energy. Therefore, 
only the electrostatic contribution to the line tension is 
analyzed. The density distribution of a conductiye liquid 
close to the three-phase contact line has been calculated 
numerically in Ref. [iE| , but the line tension was not stud¬ 
ied. Recently, Dorr and Hardt [E^ studied the electric 
double layer structure close to the TPCL by solying the 
linearized Poisson-Boltzmann equation in a wedge geom¬ 
etry, without calculati^ the line tension. Eollowing the 
method used in Ref. [E^, Das and Mitra calculated the 
Maxwell stress and the contact angle of drops or bubbles 
on a charged substrate, again without taking into account 
line effects [El| . More recently. Dorr and Hardt [48| com¬ 
puted the line tension of an electrolyte in contact with a 
charged substrate by considering a wedge geometry sim¬ 
ilar to Refs. [lEI, [EO, [El|- Similarly to Refs, [i^, 
they considered only the electrostatic contribution to the 
line tension. Howeyer, their model differs from the one in 
Refs. jlEI, 113 in that it incorporates the deformation of 
the fluid-fluid interface near the TPCL relatiye to planar 
shapes. To our knowledge there are no microscopic cal¬ 
culations of line tensions in electrolyte solutions in which 
both solyent and ion contributions are taken into account 
simultaneously. 
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Here we present a microscopic calculation of the line 
tension and of the intrinsic TP CL structure for a lat¬ 
tice model of an electrolyte solution in contact with a 
charged substrate which takes into account solvent and 
ion contributions via classical density functional theory. 
The wetting phenomena of this model have already been 
studied in Ref. 0. In Sec. mi we recall this model and 
the corresponding density functional. The results for the 
line tension and the TP CL structure for both the salt- 
free solvent and the electrolyte solutions are discussed in 
Sec. imi We conclude and summarize our main results in 
Sec. IVl 


II. MODEL AND DENSITY FUNCTIONAL 
THEORY 

A. Model 


presence of ions a first-order wetting transition occurs 
to which a prewetting line is attached from which layer¬ 
ing transition lines depart towards higher temperatures 
(see Fig. 7 in Ref. 0). These latter artifacts due to the 
lattice model used here are not expected to influence the 
results below because in the following only bulk states at 
liquid-vapor coexistence below the wetting transition are 
considered. 


B. Density functional 

We denote the dimensionless lattice positions as f = 
r/a and piiv) = pi(Y)a~^ with i G {0, -h, —} denotes the 
number densities of the solvent (i = 0) and of the ±- 
ions. The equilibrium profiles po, P+, and p_ minimize 
the following grand canonical density functional: 


We study a semi-infinite lattice model for an electrolyte 
solution in contact with a charged wall. This model is 
the same as the one used in Ref. 0. It consists of 
three components: solvent (0), anions (—), and cations 
(-h). The z axis is perpendicular to the wall. The re¬ 
gion above the wall, which is the one accessible to the 
electrolyte components, is divided into a set of cells the 
centers of which form a simple cubic lattice {r} with lat¬ 
tice constant a. The volume of such a cell corresponds 
roughly to the volumina of the particles, which are as¬ 
sumed to be of similar size. The centers of the molecules 
in the top layer of the substrate form the plane z = 0. 
At closest approach the centers of the solvent molecules 
and ions are at z = a. The plane ^ = a/2 is taken to be 
the surface of the planar wall. Each cell is either empty 
or occupied by a single particle. This mimics the steric 
hard core repulsion between all particles. Particles at dif¬ 
ferent sites interact among each other via an attractive 
nearest-neighbor interaction of strength u which is taken 
to be the same for all pairs of particles. In addition, ion 
pairs interact via the Coulomb potential. 

The wall attracts particles only in the first adjacent 
layer via an interaction potential of strength which is 
the same for all species. In addition it can carry a surface 
charge density d = aea~^ which is taken to be localized 
in the plane z = aj^ and which interacts electrostatically 
with the ions; e > 0 is the elementary charge. The surface 
charge density d is assumed to be laterally uniform and 
independent of the structure of the adjacent fluid, i.e., it 
is the same for a liquid-wall and for a gas-wall interface. 
This situation is typically realized in EWOD (electrowet¬ 
ting on dielectrics) setups [H, [53| , in which the wall is 
composed of an electrode covered by a micron-sized iso¬ 
lating dielectric layer. Eor these systems the areal wall 
charge density is determined by the laterally uniform ca¬ 
pacity of the dielectric layer and not by charge-regulation 
mechanisms or by the electric double layer structure in 
the electrolyte solution. 

Eor the present model it is known [4l| that in the 
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where [3 = (ksT) ^ is the inverse thermal energy; pi 
is the chemical potential of species i; Ib = Ibo^ = 
e^/3/(47r5o) is the Bjerrum length in vacuum; and f* = 
r*/a, p*(f*) = pi{f) for all f* G and f G with 
max(|x* — x|, 1^* — ^1, |z* — z|) < 1/2, i.e., with f cor¬ 
responding to that site of the discrete cubic lattice 
being located closest to position f* in the continuous 
space The pair potential common for all species is 
re (|r — r'l) = —u for nearest neighbors (i.e., > 0 corre¬ 

sponds to attraction) and re (|r — r'|) = 0 beyond; —u^ 
is the strength of the attractive {u^ > 0) substrate po¬ 
tential acting on the first layer z = a. D = Dea“^ is the 
actual electric displacement generated by the ions and 
the surface charge density a = aea~^, satisfying GauB’s 
law with the dimensionless gradient V obtained by 
rescaling with a: 


V ■ D (r, [pl]) = - 1/2); (2) 

i 


The concept underlying this form of GauB’s law is that all 
microscopic charges besides the ionic monopoles and the 
surface charges, e.g., those due to permanent or induced 
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dipoles, are implicitly accounted for in terms of the rela¬ 
tive permittivity 5(po(r*)). Here the relative permittivity 
is assumed to be dominated by the solvent properties, as 
it is the case for polar solvents such as water, so that it 
depends only on the solvent density po(f) but not on the 
ion densities p±(f). 

The description in Eq. m does not account for the 
structure of a hydration shell, neither in the bulk nor at 
interfaces. For actual systems there is a strong depen¬ 
dence of, e.g., the value of the interfacial tension of a 
liquid-vapor interface on atomistic details of ion hydra¬ 
tion [HI; but it is not the aim of the present study to 
model such details. Moreover, here the solubility of ions 
is accounted for merely effectively via the ion-solvent in¬ 
teraction n;(|r —r'l), which, for the sake of simplicity, 
is the same between all particle species. More realistic 
descriptions could be used instead, e.g., in terms of the 
Born energy [HI , but this is not done here for reasons of 
simplicity. 

The bulk phase diagram, i.e., the solvent and the ±- 
ion densities in the liquid {{pi,i}) and in the gas phase 
of the solution, has already been determined in 
Ref. [^. The bulk equilibrium densities are calculated 
by minimizing the bulk grand canonical potential 

= Po(lnpo -mS) + 

+ {l-po- 21) In {1-po- 21) - T(po + 2I)\ (3) 

where / :=/)+ = p_ (due to local charge neutrality) is the 
so-called ionic strength for monovalent ions; /ig = 

/ij = yd(/i+ + /i-), T* = is the reduced temperature, 
and V = Va^ is the volume of the fluid. The last term in 
Eq. ([T]) vanishes because in the bulk D = 0 due to Eq. (| 2 ]) . 
For / = 0, the reduced critical temperature is T*(/ = 
0 ) = 0.5 and the critical number density of the solvent 
is po,c{I = 0) = 0.5. For / 7 ^ 0, the reduced critical 
temperature T* is independent of / whereas po,c{I) = 

0.5-2/ nn. 

At two-phase coexistence for temperatures below the 
critical point, T* < T*, the bulk densities {pi,i} in the 
liquid and {pi^g} in the gas phase are fully specified by 
the four values po,/, lu Po,g^ /^, which have to fulfill 
the three coexistence conditions (see Eq. (26) of Ref. [4l| i 

Mo[Ka,T*]=/io[K,},T*], 

PI[{p^A.T^]=PI[{p^^g}.T% 

(4) 

where p = —Q/V is the pressure. Hence, in addition to 
Eq. (|4|), one of the four values po,z, Ih po,g^ Ig can be fixed 
arbitrarily; in the following we choose the ionic strength 
in the liquid, Ii, which, is simply called the ionic strength 
/ in Sec. [ml 

With the bulk properties fixed and with a given expres¬ 
sion for the relative permittivity 5 (po(r *))5 fbe functional 
in Eq. o can be used to study the wetting behavior of 



FIG. 1. Generic shape z = i{x) of the liquid-gas interface 
(full line) near the three-phase contact line (TPGL) on a ho¬ 
mogeneous planar surface; io is the equilibrium thickness of 
the microscopicly thin liquid-like film for T below the wet¬ 
ting transition temperature Tw and at two-phase coexistence. 
The dashed lines represent the asymptotes of ^ = i{x) for 
X ^ 00 and X ^ — 00 . The position x = 0 is defined as 
the point at which the asymptotes intersect. The asymptote 
z — i{x ^ 00 ) intersects the substrate with the contact angle 
P. The density distribution is translationally invariant along 
the y direction perpendicular to the x-z-plane. 

electrolyte solutions as it has been done in Ref. [4l| . 
If the substrate potential depends only on the direction 
z perpendicular to the wall and if the boundary condi¬ 
tion for z ^ 00 is laterally homogeneous, the equilib¬ 
rium density profiles depend on z only but not on x or 
y. Accordingly the grand canonical functional in Eq. (HD 
decomposes into the bulk contribution given by Eq. m 
and into a surface contribution proportional to the area 
A = Aa^ of the substrate: 

= VQ.b[pi]+ An,[{pi{z)}]. ( 5 ) 

The substrate-liquid surface tension and the 

substrate-gas surface tension ^s,g are given by 

Is,I = min for pi{z oo) = pi^i, (6) 

and 

^s,g= min Os[{pi(z)}], ioT pi{z ^ oo) = pi^g, (7) 
{Pi{z)} 

respectively. The equilibrium density profiles which min¬ 
imize Eq. (ED correspond to the substrate-liquid inter¬ 
face, while the ones which minimize Eq. © correspond to 
the substrate-gas interface, which consists of the emerg¬ 
ing substrate-liquid and liquid-gas interfaces separated 
by a liquid-like layer of thickness -^o(^) = CLio{T). At 
two-phase coexistence, i.e., for //* = /i*co: both density 
profiles described above are equilibrium density distribu¬ 
tions. 

Imposing these two distinct boundary conditions at 
X = Too for z ^ oo, i.e., pi{x oo^z ^ oo) = piy 
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and pi{x —oc,z ^ oc) = pi^g^ the minimization 

of Eq. ([1]) leads to an equilibrium density distribution 
pi{x^ z) which interpolates smoothly between a substrate- 
gas interface at x ^ — cxo and a substrate-liquid inter¬ 
face at X ^ +00. A specific definition of the local posi¬ 
tion of the liquid-gas interface renders a curve z = i{x) 
(see, c.f., Eq. ([IS])) such that i{x — oo) = io{T) and 


£{x ^ oo) = io{T)-\-xtdiiii^, where'd is the contact angle 
(see Eig. [1]). This arrangement leads to the formation of 
a straight TP CL independent of y where the liquid-gas, 
the substrate-gas, and the substrate-liquid interfaces in¬ 
tersect. 

Eor pi{f) = pi{x,z), the density functional in Eq. ([T]) 
can be written as 


l3^[{Pi{x,z)}] ^ I -Y^p^(^x,z)^ In (lpj{x, z)j 

x=—nxz=l \ i i j 

~ ^^Piix,z) {pj{x+l,z) + pj{x-l,z) + Pj{x,Z + l) 

+pj{x,z-l) + 2pj{x,z)) - pu^'^Pi{x,z)5zp - /3'^piPi{x,z)^ 

pLx/2 ni 

- / 

J-L^/2 Jl/ 




£{Poix,z)) 


where with lateral system size L^, L = L/a is Clausius-Mossotti expression 

the contact line length in the invariant y direction, and 
the fluid volume is E = L^L^L; = 0, +, —. 

GauB’s law (Eq. |2|) can be written as 


i-^y,{x*rz*y 


(8) 


(11) 


V-^{x\z\[pl])=^qiP*{x*X). (9) 

with the boundary conditions 

which follow from the overall charge neutrality. 

The relative permittivity 5(^*) is taken to depend 
locally on the solvent density Po(^*5^*) through the 


where a is an effective polarizability of the solvent 
molecules. In the following its value is chosen such that 
5 = 60 for po = 1; this choice corresponds to a mean value 
for liquid water along the liquid-gas coexistence curve. 
It is well-known that the Clausius-Mossotti relation be¬ 
tween the relative permittivity 5 and the polarizability a 
holds only for dilute gases. However, Eq. m is merely 
used as a simple functional form in order to obtain the 
dependence on pQ of the relative permittivity with a be¬ 
ing a fitting parameter which is adapted to interpolate 
between ^(po) ^ Po ^ ^ large value (here 

^(Po) 60) for po ^ 1. 

The Euler-Lagrange equations, which follow from the 
minimization of Eq. ^ analogously to the procedure 
presented in Subsec. IIC in Ref. [4l|, are given by 


lnpi(^>^) - Pi - /3 mw<^2,i - In (l - ^pj{x,z)^ 

3 

~ ^ E ^) + Pjix - 1, ^) + Pj{x, z+l)Jt pj{x, z - 1)) 


■ft 


l'Z+l/2 fZ+1/2 (r>(f* z* 

/ dz*4>{x*,z*)-2TrlB dz* ^^ ) y e' {pl{x*,z*))5i,o = 0, (12) 

Jz-1/2 Jz- 


lz-i/2 {e{pl{x*, z*))y 


with i,j = 0,-h, —, where qie is the electric charge of component i, T* = is the reduced temperature and 
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111 = P/ii. 0(x*,z*) = y5e0(x*,2;*) is the dimensionless 
electrostatic potential which fulfills 

D(r,r) = -^v</.(r,r). (13) 

At the wall the convention pj{x^ 0) = 0 is used. 

For given chemical potentials pi^co at coexistence, the 
coupled equations in Eq. (Ha are solved for {pi{x^z)} 
numerically by applying a Picard iteration scheme. The 
electrostatic potential is calculated for each it¬ 

eration step by solving Poisson’s equation 

V • {e{pl{x*,r))V4>{x\-z*)) = -A^Ib Y. q^P*{x* ,-z*), 

(14) 

which is a nonlinear integro-differential equation for 
after eliminating ^*) by means of Eq. ([12]). 

C. Line tension calculation 

The line tension r is calculated from the equilibrium 
density profiles {pi{x,z)} using the following definition 
for r: 

+' ' ' 5 

a=q,l 

(15) 

where Va is the volume of phase a with a G {g,l} and Qa 
is the bulk free energy density of this phase; js,g^ ls,h arid 
are the interfacial tensions and and Ai^g 

the corresponding inter facial areas of the substrate-gas, 
substrate-liquid, and liquid-gas interfaces, respectively. 
L is the length of the three-phase contact line, r is the line 
tension and • • • denotes subleading terms which vanish 
for macroscopically long contact lines L ^ 00 . 

The plane z = 0 is chosen as the substrate-fluid di¬ 
viding interface. In Ref. [T8| it has been proposed that 
in order to determine the line tension r unambiguously 
from microscopic calculations in a finite box, its bound¬ 
aries have to be chosen such that the interfaces are cut 
perpendicularly and that its edges are placed inside the 
homogeneous regions of the system. Here, in order to 
calculate the line tension, the integration box proposed 
in Ref. has been used (see Eig. 7 in Ref. |l8| and, 
c.f., Eig.[T2|). However, in a lattice model this type of box 
introduces technical difficulties for the integration proce¬ 
dure which lead to numerical errors (see Appendix for 
more details). Therefore, in order to verify the consis¬ 
tency of the results, the line tension has been calculated 
for various sizes of the integration box as described in 
Appendix 


D. Choice of parameters 

The values of the parameters used here are the same 
as the ones used in Ref. [4l|. The lattice constant a is 



FIG. 2. Local positions of the liquid-gas interface for two 
temperatures with u^/u — 0.69 for the pure solvent. The 
system undergoes a second-order wetting transition at ^ 
0.95 X T* (see Fig. 2(b) in Ref. 0). The positions of the 
liquid-gas interfaces z = l{x) (full lines) have been calculated 
from the density profiles pq{x^z) according to Eq. (1161) . For 
both X ^ 00 and x ^ — 00 , the curve z = i(x) approaches the 
asymptotes (dashed lines) from above. The position x = 0 is 
defined as the point at which the asymptotes intersect. 

chosen to be equal to 4 A, so that the maximal density 
1/a^ lies between the number densities for liquid water 
at the triple point and at the critical point. Accord¬ 
ingly, the choice Ib = 100 corresponds to T ^ 417 K. 
This temperature lies between the triple point tempera¬ 
ture of 273 K and the critical point temperature of 647 K 
for water. In our units 1 mM = 10“^ mol/liter corre¬ 
sponds to pi = piO^ = 3.9 X 10“^. The values for the re¬ 
duced surface charge density a are in the range between 
0 and 10“^. For a = 4 A the latter value corresponds to 
l/rC/cm^, which can be achieved for an EWOD setup 
(see Sec. Ill Ap by applying the moderate voltage of 60 V 
across a 100 nm thick isolating dielectric layer with a typi¬ 
cal dielectric constant eiayer = 2. In these units (Sra = 0.1 
corresponds to r ^ 1.4 x 10“^^ N. 


III. STRUCTURE OF THE THREE-PHASE 
CONTACT LINE AND LINE TENSION 

A. Line tension of the pure solvent 

First, we consider the case 7 = 0. As explained in 
Refs. 0, [53, Ell, in this case the ratio u^/u = 3T*/3i4w 
controls the wetting and drying transitions. For u^/u > 
1 the substrate is so strong that it is already wet at 
T* = 0; within the range 0.5 < u^/u < 1 there is a 
wetting transition at > 0; and within the parameter 
range 0 < u^/u < 0.5 a drying transition occurs. Here, 
the liquid-gas interfaces near the TP CL and the line ten¬ 
sion are studied for the specific choice u^/u = 0.69, for 
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FIG. 3. Dependence of the line tension r on the contact angle 
for the same pure system as in Fig. [21 The numerical results 
for the line tension (•) are consistent with the predictions of 
Ref. [H for systems exhibiting second-order wetting transi¬ 
tions with short-ranged interactions, i.e., r is negative and 
for -d ^ 0 it vanishes as r ~ —-d. The uncertainty of the line 
tension values corresponds approximately to the size of the 
symbols. The dashed line is a fit. For details concerning the 
calculation of the line tension see Appendix lAl 

which the system undergoes a second-order wetting tran¬ 
sition (see Fig. 2(b) in Ref. [4l|) at 0.95T*. We 

note that second-order wetting transitions in a pure sol¬ 
vent with short-ranged interactions is not very realistic as 
most wetting transitions either are of first order (due to 
weak van der Waals interactions) or comprise a first-order 
thin-thick transition followed by a second-order wetting 
transition (due to strong van der Waals interactions) [59| . 
However, the order of wetting transitions of the pure sol¬ 
vent is not important here; instead we intend to exploit 
the technical advantages offered by short-ranged interac¬ 
tions for the present study. 

Figure [2] shows the temperature dependence of the 
shape i{x) of the local liquid-gas interface position de¬ 
fined as 


FIG. 4. Shapes of liquid-gas interfaces for an electrolyte so¬ 
lution, which exhibits a first-order wetting transition, for var¬ 
ious surface charge densities a = crea“^, fixed temperature 
T* = 0.8 X T*, u^/u = 0.69, and fixed / = 3.9 x 10“^ 
(/ = ImM). Note that for the electrolyte solution the wet¬ 
ting temperature Tw depends on the surface charge den¬ 
sity (see Fig. 5 in Ref. [4ll]b The shapes £{x) (full lines) 
have been obtained from the density profiles po{x,z) using 
Eq. ([T6l) . For first-order wetting, the local interface position 
z = i(x) approach their asymptotes (dashed lines) from above 
for x ^ — oo and from below for x ^ oo. The position x = 0 
is defined as the point at which the pair of asymptotes inter¬ 
sect. 


this case, the line tension r is negative and vanishes as 
r ^ — 1 ^. The order of magnitude of /3ra ^0.1, which 
corresponds to r ^ 1.4 x 10“^^ N, is comparable also 
with values obtained from other theoretical approaches 
for one-component^harge-free fluids and from com¬ 
puter simulations [I, [9| as well as with experimental re¬ 
sults [3, [11-0 • 


B. Line tension of an electrolyte solution 


Lz 

'^{po{x,z) - po,g) 


In the case of second-order wetting transitions, the curve 
z = i{x) approaches the asymptotes for x ^ oo and 
X —oo from above (Fig. [2]). This result is in qualitative 
agreement with previous ones also obtained in the pres- 
ence of second-order wetting transitions 

The line tension as a function of the contact angle 
1 ^ is presented in Fig. |3l The contact angle has been 
changed by varying the temperature T*. The results 
for the line tension are compatible with the prediction 
of the interface displacement model (IDM) [Tl| for a sys¬ 
tem with short-ranged interactions approaching a second- 
order wetting transition at two-phase coexistence. In 


In this section we study the influence of the ionic 
strength / = Ia~^ and of the surface charge density 
a = aea~^ on the TPCL and the line tension. As dis¬ 
cussed in Refs. [4l|, |42|, within the chosen lattice model 
for an electrolyte solution, if cr 7^ 0 and / 7^ 0 the sys¬ 
tem undergoes a first-order wetting transition, irrespec¬ 
tive of the order of the wetting transition of the pure 
solvent. In this case, the wetting transition tempera¬ 
ture decreases with increasing surface charge density 
a of the substrate for fixed ionic strength / or with de¬ 
creasing ionic strength / for fixed surface charge den¬ 
sity a. Therefore, there are three different routes to 
vary the contact angle: (i) changing the reduced tem¬ 
perature T* and keeping the surface charge density a 
and the ionic strength I fixed; (ii) changing the surface 
charge density of the substrate a and keeping the temper- 
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ature T* and the ionic strength / fixed; and (iii) chang¬ 
ing the ionic strength I and keeping the temperature T* 
and the surface charge density a fixed. Here we con¬ 
sider the routes (i) and (ii) for two values of the ionic 
strength: / = 3.9 x 10“^ (/ = ImM) and / = 3.9 x 10“"^ 
(/ = lOmM) with Uy^/u = 0.69. 

Figure m shows the shape of the liquid-gas interface as 
obtained from Eq. (p!6]) for fixed temperature T* = 0.8T*, 
fixed ionic strength / = 3.9 x 10“^ (/ = ImM), and 
for three different values of the surface charge density a 
(route (ii)). If the wetting transition is first order, the 
local interface profile z = i{x) approaches its asymptote 
from below for x ^ oo and from above for x ^ — oc. For 
large contact angles, i.e., for small values of a (which is in 
line with the corresponding statement at the beginning 
of the previous paragraph), in Fig. SI z = £{x) follows its 
asymptotes closely. The deviation from the asymptotes 
increases for decreasing contact angles. The behavior of 
the shape of the liquid-gas interface is similar for the case 
in which the contact angle is changed using route (i). 
These results for the shape of the interface are in line 
with those of Refs. [H, [III, [29|, [6^ for first-order wetting 
in charge-free fluids. 

Figure [5] shows the line tension for the case in which the 
contact angle is changed using route (i) for two distinct 
values of the ionic strength I and for a constant surface 
charge density a = 1 x 10“^ (d = 0.1/iC/cm^). Accord¬ 
ing to Fig.infa), below the wetting transition the line ten¬ 
sion r(T*,cr,/) < 0 is a monotonically increasing func¬ 
tion of the temperature T* < T^(cr,/). Consequently, 
since T^(cr, /), and thus the deviation T^(cr, /) — T* from 
the wetting temperature, increases upon increasing / 
for fixed a (see Fig. 5 in Ref. [IJ), the line tension 
< 0 decreases upon increasing I for fixed T* 
and cr. Moreover, as discussed in Fig. [6fa) below, the 
line tension ^ a^ I) < 0 is a monotonically increasing 

function of the surface charge density a > 0 for fixed 
T* and /. Therefore, the line tension r(T*,cr,/) of the 
pure, salt-free (/ = 0) solvent in contact with a neutral 
(cr = 0) wall, which is also shown in Fig. O can be larger 
or smaller than the one for the cases / > 0, cr > 0. 

The line tension is negative and its strength decreases 
upon decreasing the contact angle (see Fig. [5[W), which 
is in line with the predictions of the IDM for the 
case of first-order wetting transitions for charge-free flu¬ 
ids with short-ranged interactions. The absolute value of 
the line tension is larger for the higher ionic strength 
I = 3.9 X 10“^ (/ = lOmM) at fixed temperature. 
We have not considered smaller contact angles because 
they require larger system sizes and therefore generate 
substantially higher computational costs. According to 
Ref. [m, the line tension in the case of first-order wet¬ 
ting transitions of fluids with short-ranged interactions 
are expected to change sign from negative to positive 
upon decreasing the contact angle'd and to be positive 
at the wetting transition temperature T^, i.e., for'd = 0. 
This agrees also with the results reported in Refs, 
for long-ranged forces. Our data do not allow us to con- 





FIG. 5. Line tension r as a function of temperature T* (a) and 
of the contact angle 'd (b) for cr = 1 x 10“^ and Uy^/u — 0.69. 
The two types of full symbols correspond to distinct values 
of the ionic strength I = la^ in the bulk liquid phase (• 
for I — 3.9 X 10“^ (/ = ImM) and ■ for / = 3.9 x 10“^ 
(/ = lOmM)); the uncertainty of the line tension values cor¬ 
responds approximately to the size of the symbols. For con¬ 
stant <j and T*, the strength |t| of the (negative) line tension 
T increases upon increasing the ionic strength I (see • and ■ 
in panel (a)). The open triangles A correspond to the case 
7 = 0 and cr = 0 (see Fig. [S]), i.e., they differ from the filled 
symbols • and ■ not only with respect to the ionic strength 
I but also with respect to the surface charge density cr. See 
Appendix lAl for details concerning the calculation of the line 
tension. Panel (b) shows the general trend of an increasing 
strength of the line tension, upon increasing the contact angle 
'd. Panel (c) displays the contact angle 'd as function of T* 
with I and cr fixed (route (i)). 


firm this prediction, but one can infer from the available 
data that such a change in sign is rather plausible. In this 
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FIG. 6. Line tension r as a function of the surface charge 
density a = go? je (a) and of the contact angle ^ (b) for 
T* = 0.8Tc and u^^ju — 0.69. The two types of symbols 
correspond to distinct values of the ionic strength / — la? in 
the bulk liquid phase (• for / = 3.9 x 10“^ (/ = ImM) and ■ 
for / = 3.9 X 10“^ (/ = lOmM)); the uncertainty of the line 
tension values corresponds approximately to the size of the 
symbols. T* = 0.8 x T? is below the wetting transitions, i.e., 
T* < T,I(cr, /), for both ionic strengths / and the whole range 
of surface charges g shown. The dashed lines are guides to 
the eye. The rightmost points in panel (a) correspond to the 
leftmost points in (b). See Appendix lAl for details concerning 
the calculation of the line tension. Panel (c) displays the 
contact angle -d as function of g with / and T* fixed (route 

(ii))- 


case, the asymptotic behavior of r for ^ 0 predicted 
in Ref. [HI is given by r ^ Tw + ci'i^ln'?^ + C2'?^ + 

Figure [6] shows the line tension for the case that the 
contact angle is varied by using route (ii) for two val¬ 
ues of the ionic strength 1 and for T* = 0.8T*. The 


FIG. 7. Gomparison of cos('d(cr)) for the systems discussed in 
Fig. [6] (solid lines) with the asymptotic expression given by 
Eq. (pTI) (dashed lines) derived in Refs, |4^ for tempera¬ 
ture T* = 0.8Tc , surface charge density 0 < a < 0.01 (0 < 
d < 1 /xG/cm^), and ionic strength / G {3.9x 10“^, 3.9x 10“^} 
(/ G {ImM, 10mM}). The asymptotic expressions apply up 
to surface charges g for which the dimensionless parameter 
X = \g\\J^'kIb![^^1^ attains unity, which is marked by the 
vertical arrows. 


line tension < 0 is a monotonically increas¬ 

ing function of the surface charge density cr > 0, and, 
as already discussed above in connection with Fig. [5fa), 
it increases upon decreasing the ionic strength I. Here, 
only small surface charge values (cr = 1 x 10“^ — 2 x 10“^) 
have been considered. Accordingly, small contact angles, 
which correspond to large surface charges, have not been 
studied. The technical reason for this is that in order to 
avoid contributions from the corners of the integration 
box, these corners should be located far away from all 
interfaces such that the density profiles near the corners 
attain their bulk values (see Appendix [Al and Sec. Ill Cp . 
Achieving this for small contact angles is more difficult 
in the case of the electrolyte solution than for the pure 
solvent, mainly due to the density distributions of the 
ions. Figures [8] and [9] show density distributions of the 
solvent po{x^ z) and of the ions p±(x, ^) for cr = 1 x 10“"^ 
(d = 0.01/iC/cm^) and cr = 8 x 10“^ (d = 0.8/xC/cm^), 
respectively. Both for Fig. [8] and Fig. [9l the bulk den¬ 
sities of the ions are = / = 3.9 x 10“^ (/ = ImM). 
For cr = 1 X 10 “^, in Fig. [8] one can see that for the posi¬ 
tive ions in the liquid phase the density profile attains its 
bulk value only in a small portion of the calculation box, 
which makes it difficult to use the integration box shown 
in , c.f., Fig. [T^ and to carry out the procedure described 
in Appendix for the calculation of the line tension. 
Moreover one can see in, c.f., Fig. [131 which shows ex¬ 
amples of the dependence of the estimator 
of the line tension (see Appendix on the box size for 
two different surface charge densities, that the amplitude 
of the variations of the value of i.e., the 

uncertainty of the value of the line tension, increases if 
the surface charge density cr increases. Figure |6l^a) shows 
that for small surface charge densities (cr < 1 x 10“^) 
the value of the line tension, the uncertainty of which 
corresponds approximately to the size of the symbols, 
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is, within the precision of the method described in Ap¬ 
pendix [3 independent of the ionic strength /. However, 
as the surface charge density increases, the absolute value 
of the line tension r decreases stronger for / = 3.9 x 10“^ 
(/ = ImM) than for / = 3.9 x 10“^ (/ = lOmM). This is 
related to the fact that due to screening for I = 3.9 x 10“^ 
(/ = lOmM) a larger surface charge is needed to produce 
the same contact angle as for / = 3.9 x 10“^ (/ = ImM) 
(see Fig. El^c)). Thus upon increasing /, according to 
route (iii), the contact angle'd increases and so does the 
strength of the line tension. 

Within the approximation of a field-free gas phase the 
asymptotic behavior 


cos{'d{a)) ~ cos(?9(0)) +—\hrr 


li,g V 2eJ 


- 1 


(17) 


has been derived in Refs. USUI for the case of the di¬ 


mensionless quantity y := \a\ 


27tIb 

el 


being small (y <C 


1). Figure [3 compares the curves cos('d(cr)) for the sys¬ 
tems discussed in Fig. [6] (solid lines) with the correspond¬ 
ing asymptotic form Eq. ([T7|) (dashed lines). The asymp¬ 
totic expressions are reliable up to surface charge den¬ 
sities a for which y ^ 1, marked by the vertical ar¬ 
rows in Fig. [3 Hence Eq. (ITT)) applies to small sur¬ 
face charge densities not only in the case of the electric 
field being confined to a wedge-shaped liquid phase, as 
in Refs. [4^, |^, but also in the case of a non-vanishing 
electric field in the gas phase. 


C. Density distributions close to the three-phase 
contact line 

The microscopic structure of the electrolyte solution 
close to the TP CL is illustrated via density maps in 
Eigs.[8]and[9]for / = 3.9 x 10“^ (/ = ImM), T* = 0.8T*, 
and two values of the surface charge density: a = 
1 X 10“'^(d = 0.01/iC/cm^) (see Eig.[8]) and a = 8 x 10“^ 
(d = 0.8/rC/cm^) (see Eig. [9]) . The contact angles are 
^ 47.5° for cr = 1X10“^ and ^ 28.3° for a = 8 x 10“^. 

Apart from the difference in contact angle and from the 
different densities of anions and cations in the vicinity 
of the wall due to the difference in surface charge den¬ 
sity cr, one can infer that for larger values of the sur¬ 
face charge the density distributions of the ions differ 
significantly from their bulk values over larger distances 
from the substrate. The anion densities p-{z) close to 
the gas-wall interface in Eig. [9] are large because in the 
present study (see Sec. Ill Ap we assume a laterally uni¬ 
form surface charge density, which is not modified by 
charge regulation. Eor setups with surface charge densi¬ 
ties being determined by charge regulation, significantly 
smaller surface charge densities would occur and hence 
smaller ion densities p±{z) close to the gas-wall interface. 

Eigures [TOl and [TTl show the charge density pc{x^z) = 
p-^{{x,z) — p-{x,z), the local ionic strength I{x^z) = 


^ (p+((x, z) -1- p_((x, 2 :))), and the electrostatic potential 
(j){x^z) = l3e(j){x^z) for the same set of parameters as in 
Eigs. [8] and [9l respectively. Eor a small surface charge 
density a (see Eig. uni the charge density pc{x^z) has a 
region in the gas close to the liquid-gas interface where 
pc{x^z) is less negative than pc{—(yz^z). If one takes a 
path parallel to the surface at small z from the gas side, 
the charge density pc{x^ z) is quasi constant in the gas 
phase far away from the liquid-gas interface, increases 
upon approaching the liquid-gas interface from the gas 
side, drops to a rather low value on the liquid side of the 
liquid-gas interface, and ultimately increases towards a 
constant value in the liquid phase. This charge separa¬ 
tion in the vicinity of the liquid-gas interface and of the 
TP CL is caused by the variation of the local permittiv¬ 
ity of the solvent e{po{x^ z)) which is higher in the liquid 
phase (see Eq. [11]). On the other hand, the structure 
of the local ionic strength distribution /(x,z), which for 
constant ^ interpolates from the value in the gas phase 
to the value in the liquid phase, is almost independent of 
z within each phase. The electrostatic potential (j){x^z)^ 
which is related to the charge density through Poisson’s 
equation (Eq. ([TT|) b does not follow the liquid-gas in¬ 
terface in that the equipotential lines bend away from 
it. Moreover, there is an electrostatic potential differ¬ 
ence between the liquid and the gas phase in the vicinity 
of the TP CL. Eor a large surface charge density a (see 
Eig. [TT]) . the high charge density pc in the vicinity of the 
substrate on the gas side screens the surface charge of the 
substrate within a few layers. This is in contrast to the 
case of small surface charge for which the charge density 
pc{x^ z) approaches its vanishing bulk value more slowly 
(compare Eig. [T0|). This different behavior is due to the 
nonlinear character of Poisson’s equation (see Eq. (HI); 
for small values of the surface charge density a its solu¬ 
tion is close to the solution of the linearized equation in 
which the number densities of the ions decay exponen¬ 
tially to their bulk values on the scale of the Debye length 
K. of the bulk phase. In contrast, for large surface charge 
density cr both the density distributions of the ions and 
the electrostatic potential (j) deviate significantly from the 
linear solution in the vicinity of the substrate, and the 
exponential decay is only valid far away from it. Eor 
this large value of cr, the aforementioned nonmonotonic 
variation of pc{x^ z) in the vicinity of the liquid-gas inter¬ 
face from the gas side is not observed. However, pc{x^ z) 
becomes more negative in the vicinity of the liquid-gas 
interface from the liquid side. This qualitative difference 
as a function of cr in the behavior of the charge density in 
the vicinity of the TP CL results in a different behavior 
of the electrostatic potential. Eor large cr the difference 
of the values of the electrostatic potential in the liquid 
and in the gas phase is not as pronounced as for smaller 
surface charge densities (see Eig. uni. Eor all cr, far away 
from the substrate the charge density pc{x^z) and the 
electrostatic potential (j){x^ z) vanish and the local ionic 
strength attains its bulk value, here / = 3.9 x 10“^. 
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FIG. 8. Density distributions of the solvent [y9o(x,z)], the cations [p+(;r,z)], and the anions [p-{x^z)] for a = 1 x 10“^ 
(d = 0.01/iC/cm^). The bulk values of the density distribution of the cations and the anions are p± — I — 3.9 x 10“^ 
(/ = ImM); pg = 0.14 and pi — 0.86 are the bulk values for the gas and the liquid, respectively. The contact angle is d 47.5°. 
The substrate is positively charged. Therefore there is a high (low) density of negative (positive) ions in its vicinity. The panels 
on the right show close-ups of the plots on the left in the vicinity of the TPCL. 



FIG. 9. Same as Fig. [8]for cr = 8 x 10 ^ (d = 0.8/iG/cm^). The contact angle is'd ^ 28.3°. Note that in this case the density 
distribution of the cations in the liquid phase needs more space in order to attain its bulk value I — 3.9 x 10“^ (/ = ImM) 
than for cr = 1 X 10“^ (d = O.Ol/xG/cm^) (see Fig. [8]). 
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FIG. 10. Charge density pc{x,z) = p-^{x,z) — p-{x,z), local ionic strength I{x,z) = 1 (p+(x,z) + p-{x,z)), and electrostatic 
potential z) = ^e^{x,z) for a = 1 x 10“^ (a = O.Ol/xC/cm^). The contact angle is ^ 47.5°. The red line in the bottom 

panels for the electrostatic potential indicates the shape i{x) of the liquid-gas interface obtained from Eq. (ITbl) . The panels on 
the right show close-ups of the plots on the left in the vicinity of the TPCL. 
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IV. CONCLUSIONS AND SUMMARY 

We have investigated the line tension and the structure 
of the three-phase contact line (Fig. [1]) of an electrolyte 
solution in contact with a charged substrate by u sing den¬ 
sity functional theory applied to a lattice model . For 
the pure, i.e., salt-free solvent, the equilibrium shape of 
the liquid-gas interface approaches its asymptotes from 
above, as expected for systems exhibiting second-order 
wetting transitions (Fig. [2j). Near the wetting transi¬ 
tion the line tension vanishes proportional to the contact 
angle (Fig. [3]) which itself goes to zero at the wetting 
transition temperature. For the electrolyte solution, the 
equilibrium shape of the liquid-gas interface approaches 
its asymptote from below as expected for systems exhibit¬ 
ing first-order wetting transitions (Fig.|4]). If the contact 
angle is changed by varying the temperature while keep¬ 
ing the surface charge fixed, the line tension becomes 
less negative as the temperature is increased (Fig. [5fa)), 
i.e., as the contact angle is decreased. For fixed tem¬ 
perature, the line tension is more negative for the larger 
ionic strength (Fig.[5ja)). If the contact angle is changed 
by varying the surface charge density at fixed tempera¬ 
ture, the line tension becomes less negative as the surface 
charge is increased (Fig. [6l^a)). For small surface charges 
this decrease of the strength of the line tension depends 
only weakly on the ionic strength (Fig. [6l^a)). However, 
for larger surface charges the decrease of the strength of 
the line tension is steeper for the smaller ionic strength 
(Fig.E^a)). We have also calculated the intrinsic equilib¬ 
rium structure of the three-phase contact line for various 
charge densities. For large surface charge densities, non¬ 
linear effects of the Poisson-Boltzmann theory dominate. 
This results in distributions of the ions and of the electro¬ 
static potential which differ from those for small surface 
charge densities (Figs. [H [H [lOl and fTT]) . 

On the one hand, technically the lattice model facili¬ 
tates the reliable determination of these structures and 
properties. On the other hand, using a lattice model 
causes a difficulty for calculating the line tension, be¬ 
cause within this model the liquid-gas surface tension 
depends on the orientation of the interfacial plane rela¬ 
tive to the underlying lattice. Accordingly, this aspect 
of our study should be regarded as a first step towards 
the microscopic calculation of line tensions in electrolyte 
solutions and should be compared with not yet available 
results from continuum models for electrolytes. More¬ 
over, for technical reasons the asymptotic behavior of 
the line tension upon approaching the wetting transition 
and the influence of a large surface charge densities of 
the substrate on the line tension could not be addressed 
within the present approach; they deserve to be analyzed 
in the future within continuum models. 


Appendix A: Line tension calculation within the 
lattice model 

For the line tension calculation, computational boxes 
B have been used which cut perpendicularly through all 
interfaces and which are bounded by the substrate-fluid 
interfaces bein g lo cated at 2; = 0 (Fig. [T^ aH. As dis¬ 
cussed in Ref. [l8|, this type of boxes ensures that, for 
sufficiently large B and within continuum models, no ar¬ 
tificial contributions to the grand canonical free energy 
Q{B) appear, which are due to the edges of B or due to 
inhomogeneities caused by the boundaries of B. Accord¬ 
ing to Eq. the grand canonical free energy fl{B) of 
B per length L of the straight three-phase contact line T 
(see Fig. [El is given by 

7s,gLs,g{B) (AI) 

where U5 = = U/ = —p is the density of the bulk 

grand potential, i.e., the negative pressure, given by 
Eq. (|3]) evaluated at the equilibrium densities, A{B) is 
the cross-sectional area of S, such that V{B) = A{B)L is 
the volume of the fluid inside S, Li^g{B) = zp{B)l sin('d) 
is the length of the intersection of the liquid-gas interface 
inside B with the x-2;-plane (see the thick magenta line 
TP in Eig. [12]), and Ls^i{B) and Ls^g{B) are the linear 
extensions of the substrate-liquid and the substrate-gas 
interface in the x-direction, respectively. The substrate- 
liquid surface tension js,i and the substrate-gas surface 
tension js,g ia Eq. (EH) do not depend on B and they 
can be inferred from the substrate being in contact with 
the bulk liquid and bulk gas, respectively. Note that the 
quantities A{B), Ls,i{B), ^s,g, and Ls,g{B) depend 
on the choice of the convention for the substrate-fluid in¬ 
terface position (here z = 0, see Eig. fl^a)). so that the 
line tension r in Eq. EH depends on this choice of the 
convention, too. 

A difference between continuum and lattice models 
arises with respect to Ji,g{B) in Eq. EH; Within con¬ 
tinuum models, Ji,g{B) = is independent of B and 

it coincides with the liquid-gas interfacial tension ")\^g ^ 
whereas within lattice models ji^g{B) varies with B since 
the tilted free liquid-gas free interface (see the thick ma¬ 
genta line TP in Eig. [12]), which is inclined by the contact 
angle'd with respect to the substrate, in general does not 
match the underlying lattice grid. 

In order to estimate the line tension r in Eq. (EH, the 
contribution ^i^g{B)Li^g{B) in Eq. EH is written in the 
form 

ji,,{B)p,{B) = (A2) 

with'd independent of B. Being maximally ignorant of 
the relative position of the liquid-gas interface with re¬ 
spect to the lattice grid, the probabilities of finding pos¬ 
itive or negative deviations Si^g{B) are equal such that 
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FIG. 12. (a) Underlying geometry for calculations of the line tension. The Euler-Lagrange equations in Eq. m are solved 
via an iterative algorithm in a rectangular box which includes as a subset the computational box B (dashed lines) which is 
used to determine the line tension. B encloses the three-phase contact line T and its boundaries cut perpendicularly through 
all interfaces. It is characterized by the lengths Ls,g and Lg^i of the substrate-gas and the substrate-liquid interface at the 
plane z = 0, respectively, by the width Lz in ^-direction, and by the ^-coordinate zp of the point P, where B intersects the 
liquid-gas interface (thick inclined magenta line, forming the contact angle with the substrate surface). The convention of the 
substrate-fluid interfaces being located at z = 0 affects not only the geometrical quantities and A{B), but also the 

definition of the substrate-gas and the substrate-liquid surface tensions 7 s,^ and 7s,z, respectively. The system is translationally 
invariant in the y-direction. Note that the boundary cutting through the liquid-gas interface crosses some of the cells, leaving 
only a fraction of each of them inside B; it is this fraction with which these cells contribute to the total grand canonical free 
energy, (b) The line tension is determined via the estimator T(B^^\B^‘^^) in Eq. (ESI . which takes two different calculational 
boxes B^^^ and B^‘^^ as its arguments (see also Eig. [Ill). 


the expectation value {6i^g{B)) vanishes. Consequently, 
according to Eq. (EH), the quantity 


T{B) 


fl{B) 

L 


- nbA{B) - 


sin('z^) 


- ls,lLs,l{B) - ls,gLs,g{B) 


(A3) 


is expected to vary, within the set of computational boxes 
B of the type specified above, as a function of B around 
the line tension r according to 


T{B) = T + 5i,g{B). (A4) 


In principle, Eq. (|A4p facilitates to determine the line 
tension r as the S-independent “background” contribu¬ 
tion to T{B). However, T{B) depends sensitively on the 
value 1 ^ of the contact angle, which turns out to be dif¬ 
ficult to track with the necessary numerical precision. A 
possible approach to determine the line tension r without 
precise knowledge of the contact angle consists of the 
following: Consider two computational boxes and 
B^‘^^ with zp{B^^^) =: 2:1 and zp(B^‘^^) =: Z 2 . The contri¬ 
butions ^ 1 / sin('z^) from Eq. (jA3|) cancel in the combi¬ 
nation ziT{B^‘^^) — Z 2 T{B^^^) so that instead of Eq. (jA4p 


one can use the expression 

Zl - Z2 

Zl - Z2 

in order to infer the line tension r as that contribution 
to T(B^^\B^‘^^), which is independent of B^^^ and B^‘^\ 
We have calculated the expression T(B^^\B^‘^^) in 
Eq. (jA5p by fixing the intersection of box B^^^ with the 
liquid-gas interface at wall distances 2:1 G [30a, 40a] for 
the pure solvent and at wall distances zi G [70a, 80a] 
for the electrolyte solution. The size of box has 
been varied accordingly such that A2: := 2:1 — 2:2 G 
{a, 2a, • • • , 15a}. This procedure has been repeated for 
all integers zi in the corresponding intervals for the pure 
solvent and for the electrolyte solution. The values of Lg^i 
and Ls^g (see Fig. [HI) are determined via the position at 
which the asymptote of the gas-liquid interface intersects 
the plane 2: = 0. The size of the rectangular box used to 
determine the equilibrium profiles depends on the contact 
angle i.e., for smaller contact angles a larger extension 
in the x-direction is needed. For the pure solvent, as the 
smaller size we have used Lx'x Ly = 300a x 60a, while as 
the bigger size 1500a x 60a has been used. For the elec¬ 
trolyte solution a fixed box size of 400a x 100a was used. 
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FIG. 13. Data for introduced in Eq. ()A5|) for an electrolyte solution with I — 3.9 x 10 ^ (/ = ImM) at 

T* = 0.8 X T* for cr = 5 X 10“^ (a = 0.05/iC/cm^) (a) and a — 1.8 x 10“^ {a — 0.18/iC/cm^) (b). The various colors 
correspond to different computational boxes in Eq. (IXsll . i.e., they correspond to different distances zi = zp(B^^^) of the 
point P, where the boundary of the box intersects the liquid-gas interface, from the substrate (see Eig. [T^ aP. The expression 
is shown as a function of — ^2 with ^2 = zp{B^‘^^) (see Eig. EHb)). The line tension r is inferred from 

these plots as that contribution to T{B^^\ B^‘^^), which is constant, i.e., independent of zi and Az (see Eq. (EH)). 


Figure [13] shows the values of calculated for 

an electrolyte solution with / = 3.9 x 10“^ (/ = ImM), 
T* = O.ST^,Uw/u = 0.69, and for two values of the sur¬ 
face charge density: a = 5 x 10“^ (d = 0.05//C/cni^) 
(Fig. fldT a)) and a = 1.8 x 10“^ (d = O.lS/iC/cm^) 
(Fig. fTdT b)). The variation of T{B^^\ with the size 
of both integration boxes, and is clearly visible. 
Nonetheless T(S*^^\S^^^) is distributed around a specific 
value r. In order to determine this value r, which here is 
called the line tension, that value of Az is chosen which 
renders the smallest variation in T{B^^\B^^^) for differ¬ 


ent ^i; r is taken to be the mean value of the smallest 
and the largest values of T{B^^\B^‘^^) for that particular 
choice A2:. We note that the amplitude of the variations 
in T{B^^\B^‘^^) increase with increasing a, i.e., with de¬ 
creasing contact angle'd. This can be inferred from the 
different scales on axes of ordinates in Figs. [T3l a) and 
(b). The corresponding behavior is similar for the other 
values of the surface charge density considered here. If 
the surface charge density is fixed and the contact angle 
d is varied by changing the temperature, the amplitudes 
of the variations in T{B^^\B^‘^^) increase upon increas¬ 
ing the temperature T*, i.e., upon decreasing the contact 
angle d. 
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